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Abstract 

We investigate collective multipole excitations for closed shell nuclei from ^^O to ^"*^Pb using correlated 
realistic nucleon-nucleon interactions in the framework of the random phase approximation (RPA). The 
dominant short-range central and tensor correlations are treated explicitly within the Unitary Correlation 
Operator Method (UCOM), which provides a phase-shift equivalent correlated interaction Vucom adapted 
to simple uncorrelated Hilbert spaces. The same unitary transformation that defines the correlated interac- 
tion is used to derive correlated transition operators. Using Vucom we solve the Hartree-Fock problem and 
employ the single-particle states as starting point for the RPA. By construction, the UCOM-RPA is fully 
self-consistent, i.e. the same correlated nucleon-nucleon interaction is used in calculations of the HF ground 
state and in the residual RPA interaction. Consequently, the spurious state associated with the center-of- 
mass motion is properly removed and the sum-rules are exhausted within ±3%. The UCOM-RPA scheme 
results in a collective character of giant monopole, dipole, and quadrupole resonances in closed-shell nu- 
clei across the nuclear chart. For the isoscalar giant monopole resonance, the resonance energies are in 
agreement with experiment hinting at a reasonable compressibility. However, in the 1~ and 2^ channels the 
resonance energies are overestimated due to missing long-range correlations and three-body contributions. 
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I. INTRODUCTION 



One of great challenges for nuclear theory is the description of ground state properties and exci- 
tation phenomena in finite nuclei, based on realistic nucleon-nucleon (NN) interactions. A variety 
of highly accurate realistic NN potentials are presently available, e.g. Nijmegen U], CD Bonn lO], 
and Argonne VI 8 |3]. Recently, realistic NN potentials have been constructed in the framework 
of chiral perturbation theory, assuming a chiral symmetry breaking scale chosen as to maintain 
pions and nucleons as relevant degrees of freedom |4]. However, a quantitative description of 
nuclear structure properties necessitates the inclusion of additional ingredients besides the two- 
nucleon interaction. In particular, it has been noted that relativistic corrections and three-nucleon 
(3N) interactions in light nuclear systems can give sizable contributions Isl. Most of the available 
3N interactions arephenomenological, i.e., their parameters are adjusted to experimental data in 
finite nuclei BaLZHSD- So far, only the chiral approaches offer a consistent derivation of two- and 
three-nucleon forces. 

These advances strengthend the interest in ab initio nuclear structure calculations based on 
realistic potentials. Recent studies include different theoretical approaches for the description of 
ground state properties and low-lying excitation spectra in finite nuclei, e.g. Green's function 
Monte Carlo |9], no-core shell model liol, and coupled cluster method lllll . Because of huge 
computational requirements, most of these methods are limited to light nuclei. In order to access 
heavier nuclei one has to resort to approximate solutions of the nuclear many-body problem. A 
particularly appealing method is the self-consistent Hartree-Fock (HF) scheme, which has been 
very successful in the description of various nuclear ground state properties. However, those cal- 
culations were usually based on simple phenomenological interactions such as Gogny or Skyrme 
interactions 1 12 , ist ^ whose microscopic foundation is not well established, since their parameters 
are adjusted to the bulk nuclear properties. The use of bare realistic NN interactions in HF dramat- 
ically fails, resulting in unbound nuclei |14]. This is a direct consequence of strong many -body 
correlations induced by the short-range repulsion and the tensor part of the realistic potentials. 
The dominant short-range correlations cannot be described by many -body states given by a simple 
Slater determinant as in the HF scheme. Therefore, for practical applications of the HF approx- 
imation, the realistic NN interaction has to be converted into an effective interaction, which is 
adapted to the many-body model space under consideration. 

This objective can be achieved within the Unitary Correlation Operator Method (UCOM), 
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which describes short- range c entral and tensor correlations explicitly by means of a unitary trans- 




formation 11151 116LI17L 118111911 . Other methods employing unitary transformations have been de 



ni 



2211 and the Lee-Suzuki trans- 



veloped, including the unitary model operator approach (2011211 

formation J^]- In contrast to these methods, the correlation operators in the UCOM approach 
are given explicitly, allowing for the derivation of a state-independent operator form of the corre- 
lated interaction and of other relevant correlated operators. Although different by its construction, 
the effective NN interaction from the UCOM method is similar to the Viow-k potential obtained 
by using renormalization group concepts ll24ll . Both approaches provide a phase-shift equivalent 
low-momentum interaction, which is appropriate for nuclear structure calculation in simple model 

nnn 

spaces II19LI25L 12611 . First applications of the correlated interaction in Hartree-Fock and many -body 



perturbation theory from He to ""Pb show the potential of this approach 



Hi. 



In the present work we study collective excitation phenomena in atomic nuclei using a cor- 
related interaction based on the Argonne VI 8 potential. A particularly convenient method to 
investigate low-amplitude excitation phenomena is the random-phase approximation (RPA). Dif- 
ferent versions of RPA and quasiparticle RPA, based on phenomenological interactions, have been 



verv successful not on 



V for the description of giant resonances and low-lying states (e.g. Refs. 



33, 
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36|]), but also in studies of exotic nuclear structure of collective 
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In 



excitations in nuclei away from the valley of yS-stability 

HQS 

the present study, correlated realistic NN interactions are employed for the first time to investigate 
collective excitations in both light and heavy closed-shell nuclear systems. This serves as a strin- 
gent test of the UCOM framework and provides direct information about the physical properties 
of the underlying correlated realistic NN interactions. 

In Seciniwe introduce the basic formalism of the random-phase approximation in the frame- 
work of the unitary correlation operator method. In Sec.|lll|we present some fundamental tests of 
validity of the scheme: the separation of the spurious center-of-mass motion and the accuracy of 
the sum rules in comparison to the standard expressions. Section|lVlpresents the application of the 
UCOM-RPA approach for the description of collective excitation phenomena, in particular giant 
monopole, dipole, and quadrupole resonances. The role of the range of the tensor correlator and 
the impact of missing long-range correlations and three-body forces is discussed. Finally, in Sec. 
IVjwe summarize our findings. 
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II. THE RANDOM-PHASE APPROXIMATION BASED ON THE UNITARY CORRELATION 
OPERATOR METHOD 



We employ correlated NN interactions constructed within the UCOM approach, for the descrip- 
tion of small-amplitude oscillations around the nuclear ground state. In a first step we solve the 
HF equations based on the two-body matrix elements of the correlated realistic NN interaction. In 
a second step, the RPA equations are formulated in the HF single-nucleon basis. Here we outline 
the basic principles of the UCOM scheme and the HF and RPA formalism. 

A. Unitary correlation operator method (UCOM) 

The Unitary Correlation Operator Method provides an eff'ective NN interaction which can be 
directly used in nuclear structure calculations ifni [isl [igj . The central idea of this approach is 
the explicit treatment of the interaction-induced short-range central and tensor correlations. These 
correlations are imprinted into an uncorrelated many-body state |*F) through a state independent 
unitary transformation defined by the unitary correlation operator C, resulting in a correlated state 

i¥) = cm. (1) 

Even for the simplest uncorrelated state, a Slater determinant, the correlated state |*F) contains 
the relevant short-range correlations. Any expansion of I^P) in a basis of Slater determinants will 
require a huge number of basis states, i.e. a large model space. Hence, the unitary transformation 
reduces the size of the model space necessary for an adequate respresentation of the full many- 
body state. 

The correlation operator C is written as a product of unitary operators Cq and C^ describing 
tensor and central correlations, respectively. Both are formulated as exponentials of a Hermitian 
generator, 

C = CnC, = exp[ - i ^ gn,;y] exp[ - i ^ g^jj] . (2) 

The two-body generators g^ and gn are constructed following the physical mechanisms by which 
the interaction induces central and tensor correlations. The short-range central correlations, caused 
by the repulsive core of the interaction, are introduced by a radial distance-dependent shift pushing 
nucleons apart from each other if they are within the range of the core. Radial shifts are generated 
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by the component of the relative momentum q = ^[Pi - P2] along the distance vector r = Xi - X2 
of two particles 

qr = ^[q-f + ^q]. (3) 

The dependence of the radial shift on the particle distance is described by a function s{r) in the 
Hermitian generator, 

gr = ^[s(.r)qr + qrs(r)] . (4) 

The application of c,- in two-body space corresponds to a norm conserving coordinate transfor- 
mation r 1-^ R-{r)^ with respect to the relative coordinate. The radial correlation function R-ir) 
and its inverse R+(r) are related to the function s{r) in the following way iIitIi . 



±1. (5) 



For a given bare potential, the central correlation functions R+{r) are determined by an energy 
minimization in the two-body system for each iS,T) channel. In the purely repulsive channel 
(S, r) = (0, 0), an explicit constraint on the range of the central correlator is used. 



dr r\R^{r) - r) = Ir^. (6) 

Throughout this work, we use the Argonne VI 8 potential with the optimal correlators constructed 
in Ref. H. We adopt the short-range central correlator with the constraint /^^ =0.1 fm"^. 

We have verified that variations around this value have negligible effect on the ground state prop- 
erties and excitation spectra in 0^, 1~, and 2^ channels, in closed-shell nuclei across the nuclear 
chart. 

Tensor correlations between two nucleons are generated by a spatial shift perpendicular to the 
radial direction. Li practice, this can be achieved by using the "orbital momentum" operator 

qn = q-7q'- = ^[Lxr-rxL] , (7) 

where L is the relative orbital angular momentum operator. Radial momentum ^q^- and orbital mo- 
mentum constitute a special decomposition of the relative momentum operator q and generate 
shifts orthogonal to each other. The dependence of the shift on the spin orientation is implemented 
in the following way in the generator gn I16I1 : 

3 

gn = ^&(r)[((ri ■ qa)(o-2 ■ r) + (o-j • r)(o-2 • q^^)] . (8) 



The two spin operators and the relative coordinate r enter in a similar manner like in the standard 
tensor operator S12, but one of the coordinate operators is replaced by the orbital momentum q^, 
which generates the transverse shift. The size and radial dependence are given by a tensor corre- 
lation function d^{r) for each of the two 5 = 1 channels. The parameters of §{r) are determined 
from an energy minimization in the two-body system lisl . with an additional restriction on the 
correlation volume which constrains the range of the tensor correlator, 



dr r^§{r) = U. (9) 



As for the central correlators we adopt the optimal correlation functions for the Argonne VI 8 
potential determined in Ref. Iisl. Applications within the no-core shell model have show that 
~ ' ~ = 0.09 fm^ leads to the best description of binding energies in and '*He H. This 
correlator set 118 8h . which we refer to as standard correlator in the following, also provides^ a good 
description of binding energies for heavier nuclei within many -body perturbation theory [19]. In 
addition to the standard correlator, we will employ other values for the constraint on the range of 
the tensor correlator, = 0.07, 0.08, and 0.09 fm^, in order to probe its relevance for the 

description of the global properties of collective excitation phenomena in atomic nuclei. 



B. Correlated operators 

Due to the unitarity of the correlation operator, matrix elements of an operator A with correlated 
many-body states are equal to those evaluated with the correlated operator A and uncorrected 
many-body states, i.e. 

<¥| A I?') = <¥| C = <*P| A (10) 

The correlated operator contains irreducible contributions to all particle numbers. Within a cluster 
expansion of the correlated operator 

A = CUC = A"! + A^^l + JP^ + • • • , (11) 

where A^''^ denotes the irreducible n-body contribution, we usually assume a two-body approxi- 
mation, i.e. three-body and higher-order terms of the expansion are neglected. In previous studies 
it has been verified that higher order contributions due to s 
neglected in the description of nuclear structure properties 
the tensor correlations. The tensor interaction is very long -ranged and thus generates long-range 



lort-range central correlations can be 
I ni . However, this is not the case for 



tensor correlations in an isolated two-body system, e.g., the deuteron. In a many-body system, 
the long-range tensor correlations between a pair of nucleons are suppressed by the presence of 
other nucleons, leading to a screening of the tensor correlations at large interparticle distances. 
In terms of the cluster expansion this screening appears through significant higher-order contribu- 
tions. In order to avoid large higher-order contributions we have restricted the range of the tensor 
correlation function (cf. Sec. Ill All . which provides an effective inclusion of the screening effect 

H). 

Starting from the uncorrected Hamiltonian for the A-body system. 




consisting of the kinetic energy operator T and a two-body potential, the formalism of the uni- 
tary correlation operator method is employed to construct the correlated Hamiltonian in two-body 
approximation 

= + jm + y[2] ^ X + VucoM, (13) 

where the one-body contribution comes only from the uncorrelated kinetic energy T^'^ = T. Two- 
body contributions arise from the correlated kinetic energy T^^^ and the correlated potential V^^^ 
which together constitute the phase- shift equivalent correlated interaction Vucom 111 811 . 

The correlated realistic NN interaction Vucom is a good starting point for a study of nuclear 
structure. However, one should keep in mind that long-range correlations and residual three-body 
forces are not yet included. One can account for long-range correlations by a suitable extension 

o 

for the many-body space as discussed in Ref. II19I1 . The problem of effective three-body interac- 
tions, which are composed of the genuine three-body force and the three-body terms of the cluster 
expansion, remains as an objective for future studies. Therefore, in the present work we focus on 
the question to which extent the correlated NN interaction alone is suflicient for the description of 
collective excitations in finite nuclei. 



C. Hartree-Fock method with correlated realistic NN-interactions 



The correlated realistic NN potential Vucom is suitable for the use in HF calculations for the 
ground state of finite nuclei [19]. We start from a Hamiltonian which consists of kinetic energy 
and the Vucom interaction derived from the Argonne VI 8 potential including the Coulomb poten- 
tial iISQIi . More details about the UCOM-HF scheme are available in Ref. lll9ll . The center-of-mass 
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contributions are subtracted on the operator level, i.e. we employ the correlated intrinsic Hamilto- 
nian in two body approximation, 



Hint - T - Tcni + VuCOM - Tint + VuCOM • (14) 

The intrinsic kinetic energy operator reads, 

Tint = T-T,,, = ^— Vq^ (15) 
AniM ^ ■' 

'<j 

where q corresponds to the relative two-body momentum operator, and we assume equal proton 
and neutron masses and thus a reduced mass // = nif^/l. Assuming spherical symmetry, the HF 
single-particle states are expanded in a basis of harmonic oscillator eigenstates, 

\vl jmm,) = C^n^'"'"'^ \nl jmm,) , (16) 

n 

with radial quantum number n, orbital angular momentum /, total angular momentum j with pro- 
jection m, and isospin projection quantum number m,. For closed shell calculations, we restrict 
^(yijmm,) independent of m. The HF equation can be written in matrix form, 

h 

which is solved self-consistently to determine the expansion coefficients and single-particle ener- 
gies. The single-nucleon Hamiltonian, 

I' ,j',m'i n',n' 

is constructed from the m-averaged antisymmetric two-body matrix elements of the correlated 
intrinsic Hamiltonian Hint <fT^ . 

jj{ljm,l'j'm',) _ 1 V"' 

nn'M' ~ (2j+ 1)(2 / + 1) 

^ ■' ■' ' m,m' (19) 

X {nl jmm,, n'l' j'm'm[\ Hint \nl jmm,, h'l' j'm'm[), 

where corresponds to the one-body density matrix. 

The harmonic oscillator basis is typically restricted to 13 major shells, which warrants complete 
convergence of the HF results. Calculations with larger basis sizes are possible but rather time- 
consuming, because of the computational effort for evaluating the two-body matrix elements of 
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the VucoM potential. The optimal value of the harmonic oscillator length ao = ylftlniNiOo is 
determined from an explicit energy minimization for different regions in the nuclide chart. 

In Figs. [Hand El we show the UCOM-HF single-nucleon spectra near the Fermi level for ^^O 
and "^"Ca, respectively. The calculations are based on the correlated Argonne VI 8 interaction, 



using the standard correlator (Sec III All . The calculated energy levels are compared with the HF 
spectra obtained with the low-momentum NN potential V/on-zt, with two standard phenomenologi- 
cal interactions in the nonrelativistic (Skyrme) and the relativistic (NL3) ll47l] framework, and 
with experimental levels which are evaluated from the binding energies For the case of ^^O, 
the UCOM-HF single-particle spectrum is similar to the one obtained with the Viow-k potential. 
These two approaches, however, result in somewhat different spectra than the phenomenological 
models and those extracted from the experiment. For HF approaches based on realistic NN in- 
teractions, the spectra appear spread too wide in energy. The results for "^"Ca are similar, but the 
UCOM-HF spectra for neutrons and protons are slightly more compressed in comparison to the 
Viow-k case. 

The origin for this behavior is in the missing long-range correlations, higher order terms of the 
cluster expansion, and genuine three-body interactions. One way to account the long-range corre- 
lations for the binding energies and radii beyond the simple mean-field approach is the framework 
of the many -body perturbation theory (MBPT). When employed in the UCOM scheme, the sec- 
ond order MBPT recovers quite a significant part of the missing binding in UCOM-HF \\^]. An 
alternative approach to account for the long-range correlations in the ground state would be the 
inclusion of the correlations due to low-lying excited states and collective excitations within an 



RPA framework 143 



D. Outline of the UCOM random-phase approximation 

In the limit of small- amplitude oscillations of the nuclear density around the ground state, col- 
lective excitation phenomena can be studied within the random-phase approximation (RPA). We 
address the question to which extent are the UCOM-HF single-nucleon basis and the residual 
interaction based on the correlated realistic NN interaction sufficient for a description of highly 
collective excitation phenomena, such as giant resonances. Since details about the derivation of 
RPA equations are available in textbooks \q\ IsoI . we review the basic principles only briefly. 
The UCOM-HF single-particle states are used for the construction of the particle-hole {ph) con- 
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figuration space for the RPA scheme. Assuming a spherical nuclear system, the coupling to good 
angular momentum is employed. The collective excited states of multipolarity S are generated by 
the quasiboson operator, 

Q;aMlO)= (20) 
where the RPA vacuum |0) is defined by the condition. 



Qv,aMlO) = 0. 



(21) 



The quasiboson operator reads, 

^y.aJVl / ; L ph ph ph ^ ' ph \^ 

ph 

where the sum runs over the ph states of the HF single-nucleon basis, and 



(22) 



ph 



J] (jpnip, jnm.l aM)(-l)^""'""a;„,^ a,-,„„, 



(23) 



ripini, 



corresponds to the ph creation operator. One of the standard approaches to derive RPA is the 
equation of motion method using the quasiboson approximation [27], with the RPA vacuum ap- 
proximated by the HF ground state, i.e. |0) ~ \HF). The resulting set of coupled equations for 
the amplitudes and I'^f ""^ and the RPA eigenvalues ojy is given by. 



ph 



( A3 B3 

g3* j^3* 



1 
-1 



^v,g3vt 

yv,a3Vt 



(24) 



The RPA matrices for the given configuration space of the single-nucleon UCOM-HF basis are 
obtained from. 



p3 

php'h' 

b3 



<//F|[[A9f,Hi„,],A3:^1|//F) 
-{HF\ [[A^f , H,,], (-l)3-^Aj-,^] \HF), 



(25) 



php'h' ^ ' i-i- ph 

where the Hamiltonian Him includes the intrinsic kinetic energy (fTSt together with the correlated 
potential Vucom, in a consistent way with the Hartree-Fock equations. In the present study, we 
assume that the RPA vacuum is rather well approximated by the HF ground state. In principle, 
however, one would need to build excitations on the RPA vacuum, and iteratively solve the ex- 
tended RPA equations Isil . A preliminary study within an extended RPA framework indicates 
that the proper implementation of the RPA vacuum causes only small corrections in the excitation 
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spectra. In particular, the centroid energies of the strength distributions systematically decrease 
by approx. 1 MeV in the isovector channel, while the excitation energies of isoscalar modes are 
affected even less and increase by less than 1 MeV. More details will be provided in a forthcoming 
publication. 



E. Transition operators 

The response for electric multipole transitions is given by the reduced transition probabil- 
ity Id], 

B\Ed, Ji ^ Jf) = ^-l^|</||Qj||0f , (26) 

where Q J corresponds to the electric multipole transition operators. The isoscalar monopole oper- 
ator is defined as 

A 

Q'^o' = Yj^iYoo{Xi), (27) 

where x, = |x,|. The multipole (3 > 0) isoscalar (T = 0) and isovector (T = 1) operators are given 
by 

A 

Qj;? = ^^ J]xf 7aM(x/) (28) 

and 

A 

QS=eJ]Tfxf Fa3vt(x,-), (29) 

i=l 

respectively. In the UCOM framework, the operators of all observables need to be transformed in a 
consistent way. Therefore, the same unitary transformation that is used for the nuclear Hamiltonian 
has to be employed for the multipole operators entering into transition matrix elements. The effect 
of using correlated transition operators will be examined in the cases of monopole and quadrupole 
transitions. Since the unitary correlation operators act on the relative coordinates, we rewrite 
the monopole and quadrupole transition operator in a two-body form with separated relative and 
center-of-mass (cm.) contributions (cf. Ref. l52lll '). 

A 

^00° = ^x?^oo(x,) 

' ^ A (30) 
= ITJTT; Z HYooihj) + 4Xj.yoo(X,,)] 
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' , A (31) 



2(A - 1) . . 

where r,y = x, - Xy are the relative and X;y = (x, + Xj)/2 the cm. coordinates of a nucleon pair. 

The correlated multipole operator is constructed in an analogous way as the correlated operators 
of the realistic NN interaction IirI . In two-body approximation we obtain 

= Q[l] + Q[2] = Q + Q[2]^ (32) 

where Q^^^ = Q corresponds to the bare one-body transition operator, and the two-body part of the 
correlated operator reads 

~ 1 ^ 

[Q^'f = E i'^' [r'/oo(f,7)] c - rfjYooihj)] (33) 

for the isoscalar monopole transition operator, and 
~ 1 ^ 

iQ2Mf' = Z i'^' [4-^2M(r,y)] c - rfjY^Mirtj)} (34) 

for the isoscalar quadrupole transition operator. Therefore, in addition to the usual transition 
matrix elements for the bare multipole operator, 

<vi 0^3^ 10) = ^Yj i^;f + (-1)^'"'^'""^ (hmm, (35) 

we also have to include contributions from the correlated two-body part. 



(36) 



phJ[J2h' \ Jl Jl ih' 

x{(-l)a+-'^x;f<(j-,jV)/ill[Qj]f'^ll(7/,^)^2) 

+ (-i)'' Y];!{Uhh')J2\\V^f\iph')Jx)], 

for each RPA eigensolution {(x)v,X^, 7"), using the standard notation ^ = V2tJ -i- 1. The resulting 
transition strength function, including both the one and two-body contribution of the correlated 
transition operator, is given by 

B'^{E2,Ji^ Jf) = Bl{io,) 

(37) 

^ — \(v\\o:im + <'viirQ'fiL2J|in\r. 

27,: 



:^|<y||Qj||0) + {vmlt\\0)\\ 
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In order to evaluate the two-body matrix elements of the transition operator we use the relative 
spin-orbit coupled states of the form \n(LS)JM TMj), with a generic radial quantum number n, 
relative orbital angular momentum L, spin 5, total angular momentum J, and isospin T. In the 
following, the projection quantum numbers M and Mj are suppressed. For the two-body parts 
of the correlated monopole (B^ and quadrupole operator we have to evaluate general matrix 
elements of the form {n(LS)JT\clcl^f(r)Ygy[(r)cnCr\n'(L'S)J'T), where Cq and c,- indicate the 
central and tensor correlation operators in two-body space. 

As for the correlated Hamiltonian it is convenient to apply the tensor correlator to the two-body 
states and the central correlator to the operator lisi . To this end we formally interchange the 
central and the tensor correlation operators using the identity CqC,- = c^cJcnCr = CrCn- For the 
transformed tensor correlator we get 

Cn = cjcnc, = exp[-ign] (38) 

with 

gn = l [(o-i • qn)(o-2- r) + {cti ■ r)(o-2- q^)]. (39) 

The tensor correlation operator acts like the identity operator on the LS -coupled two-body states 
with L = 7 O, while for L = 7 + 1 it gives 

Cn \n{J + 1, l)7r) = cosOjir) \n{J + 1, l)JT) 

(40) 

+ %m6j{x) K7+ l,l)7r) , 

where 

?Xr) = 3 ^JJ{J+l) • (41) 

Using Eq. (HUI) we can easily derive explicit expressions for the correlated two-body matrix el- 
ements {n(LS)JT\c^iicl f(r)Ygjyi(r)ci^ii \n'{L'S)J'T) of monopole and quadrupole transition op- 
erators. For L = J and U = J' the tensor correlator is inactive and the correlated matrix element 
assumes the simple form 

{n(JS )JT\ clcl mY3M(r) c^c, \n'{J'S )J'T) = 

(42) 

dxulj{r)u„A^) m {{JS)JT\ Y3^{r)\{J'S)J'T), 



where /(r) = f(R+(r)) is the transformed radial dependence of the multipole operator, and Un^i^) 
is the relative radial wave function. For L = J + I, L' = 7' + l and 5 = 1 the tensor correlator, 
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following Eq. (BUI) , transforms the bra and the ket-state of the matrix element into a superposition 
of two states. This leads to a combination of four terms for the correlated matrix element 



{n{J+l, l)7r|cJc^/(r)yaM(f)cnC.|n'(/ + l, \)J'T) = 

= ^ drM*j^i(r) Un'j'^iir) f{r) 

X [<(7+l, l)7r| Y^Mir) 1(7' + 1, 1)7T) cos Wcos?y,(r) 
+ <(7±1, l)7r| ygM(r) 1(7' ±1, 1)7T) sin?X^) sin?,,^-) 
± <(7+l, l)7r| Y^Mix) 1(7' ±1, l)7'r) cos?,(r)sin0,.(r) 
± <(7+l, l)7r| Fg3vc(r) |(7'±1, l)7'r> sin?X^) cos 0,,(r)] 



An analogous expression is obtained for L = 7 + 1 and U = 7' + 1. For L = 7 + 1 and U = J' , the 
tensor correlator affects only the bra-state and we get the simpler form 



At this point we should like to remark, that the present UCOM-RPA scheme is self-consistent in 
two respects. First, the same correlated realistic NN interaction Vucom is used in the HF equations 
that determine the single-particle basis and in the RPA equations as residual interaction. Hence, 
the effective NN interaction which determines the ground-state properties, also determines the 
small amplitude motion around the nuclear ground state. This ensures that RPA amplitudes do not 
contain spurious components associated with the center-of-mass translational motion, which will 
be shown explicitly in Sec. |nT| Second, the unitary transformation used to construct the effective 
interaction Vucom is also applied to the transition operators. Hence the effect of short-range central 
and tensor correlations is included consistently in all relevant observables. 

ni. TESTS OF THE UCOM-RPA IMPLEMENTATION 

In order to ensure that the UCOM-RPA scheme is properly implemented, and to probe its self- 
consistency, we perform several stringent tests. In particular, in studies of the multipole response 
of closed-shell nuclei across the nuclear chart, the following conditions need to be fulfilled: (a) 



(n(7+ 1, l)7r| cjc^ /(r)yaM(r) cnc, |n'(7', l)7'r) = 




(44) 



X [<(7±1, l)7r| Y^MimJ', ^V'T) cosOjir) 



+ <(7+l, l)7r| YsMir) 1(7', l)7'r) sin ■ 
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the spurious excitation corresponding to a translation of a nucleus decouples as a zero-energy 
excitation mode, (b) the transition strength should exhaust the sum rules, and (c) the excitation 
energies of giant resonances and low-lying states converge to stable solutions independent of the 
energy cut-off parameter. As will be demonstrated in this section, the present UCOM-RPA scheme 
is fully consistent with these conditions. 



A. Spurious solutions of the UCOM-RPA equations 

Whenever the generator of a continuous symmetry for a general two-body Hamiltonian does 
not commute with the original single-particle density, it produces a spurious zero-energy solution 
of the RPA equations Isoll . It corresponds to a mode which is not related to an intrinsic excitation 
of the system, but in fact to a collective motion without restoring force. In our particular case, a 
spurious dipole mode is associated with translation, i.e. the center-of-mass motion. Ideally, if the 
RPA scheme is built on the self-consistent wave functions and single-particle energies, the spurious 
excitation should decouple from the physical states at exactly zero energy. In practical calculations, 
however, due to the truncation of the ph configuration space and inconsistencies between the 
ground state and RPA equations, the spurious state is separated at energies larger than zero. In 
this case the physical states may be more or less mixed with the spurious response, leading to 
seriously overestimated strength distributions, especially in the low-energy region [33]. Violations 
of self-consistency in the Hartree-Fock based RPA may also cause a spurious enhancement of the 



isoscalar monopole mode energy for spin unsaturated systems 115311 . 

The truncation of the actual RPA configuration space may strongly aff"ect the separation of 
spurious states. The size of the ph space is determined by a cut-off parameter Eph-MAx which 
corresponds to the maximal allowed energy for ph excitations. In Fig.|3]we display the energy of 
the spurious state Ess as a function of the energy cut-off Epi^^MAx, for '^O and ^"*^Pb. The isoscalar 
dipole operator is employed, 

A 

Qj;^ = eJ]x^Y:MiXi), (45) 

which leads to both, Ihco and 3fia) ph excitations. In the same graphs, we also present the corre- 
sponding number of ph configurations Npf, which enter into the RPA equations. We employ the 
correlated Argonne VI 8 interaction with the standard correlators, as discussed in Sec. Ill Al As 
the number of ph configurations increases, the spurious state converges towards zero excitation 
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energy for both nuclei. This necessitates relatively large spaces with 80 and 330 configurations 
for '^O and ^°^Pb, respectively. The accuracy of our method is exemplified by the energies of 
the spurious states : 0.01 MeV for ^^O, and 0.05 MeV for ^°^Pb, which appear to be lower than 



in recently deve 

r 

interactions 1138 



oped fully self-consistent (Q)RPA models based on phenomenological nuclear 



Next, we employ eff'ective dipole transition operators, 

A 



= ^ Z - ^ (x')ox,) Y,m(%X (46) 



A 



^;l=eJ]<^^-^^^iY,^d, (47) 

n 

which explicitly contain the center-of-mass correction terms on the operator level II54I1 . The high 
accuracy of the UCOM-RPA method is illustrated in Fig. HI for "^^^Ca, where we compare the 
isoscalar dipole transition strength distribution for the uncorrected operator (l45t . and the oper- 
ator with the center of mass correction term The two spectra are almost identical, and the 
major difference is only in the strong transition of the spurious state at 0.005 MeV, obtained for 
the operator without a center of mass correction. When the correction term is included, the spu- 
rious state is completely removed, whereas transitions corresponding to actual dipole vibrations 
are practically unaffected. This means that the physical excited states are free of spurious con- 
taminations, owing to the full self-consistency of our method. In Fig.lHone can also observe a 
double hump structure that is characteristic for the isoscalar giant dipole resonance (ISGDR). It 
is composed of a high-energy part corresponding to a compression mode, and a low-energy part 
which might be of different nature ll55h . 



B. Sum rules of the multipole strength in the UCOM-RPA scheme 



The sum rules of various collective excitation modes is of particular interest, since their values 
represent a useful test for the RPA ll50ll . In the present study, we examine the energy- weighted 
sum rule (EWSR) 11381] for the isoscalar monopole. 



and isoscalar quadrupole excitations. 



(48) 



(49) 
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In Fig. |5l we plot the cumulative energy weighted sum of the transition strength for given 
excitation energy E < 50 MeV in the closed-shell nuclei ^^O, "^^Ca, ^°Zr, and ^"^Pb. The upper 
and lower panels display the summed energy weighted strength for the isoscalar monopole and 
quadrupole cases, respectively. One can observe that beyond 40 MeV the sums converge to their 
final values, which are in good agreement with the EWSR from Eqs. (BHl) and ( B^l) . We have 
also confirmed that UCOM-RPA essentially exhausts the EWSR with maximal discrepancies of 
±3% in other closed-shell nuclei. Accordingly, we conclude that the completeness properties and 
consistency are accurately fulfilled in the UCOM-RPA approach. 

In the following, we will average the discrete UCOM-RPA strength distributions with a 
Lorentzian function, and therefore obtain the continuous strength function. 



The width of the Lorentzian distribution is fixed to the arbitrary value Y=2 MeV. The Lorentzian 
function (l50b is defined in a way to fulfill the condition that the sum of the energy weighted 
response is equal for the discrete distribution and the continuous strength function. 



C. The role of correlated multipole transition operators 

The UCOM-RPA observables describing collective excitation phenomena are evaluated by a 
consistent application of the same unitary transformation as for the nuclear Hamiltonian. In this 
section, we test the relevance of the correlated transition operators for the multipole strength dis- 
tributions. 

In Fig. |6l the isoscalar monopole and quadrupole responses in ^^O are displayed for the two 
cases: (i) with the bare multipole operators dTTl . (EHb . (l35t . and (ii) with the correlated multipole 
operators constructed by the unitary transformation, Eq. (071 . The two transition strength distri- 
butions are essentially identical. In order to probe the relevance of the tensor correlations for the 
quadrupole response (in the monopole case, only the central correlations are active), we construct 
a new tensor correlator with a long range, constrained by ' =0.2 fm . However, even in this 
extreme case, the correlated quadrupole response resembles the one obtained with the bare oper- 
ator (Fig.|6l). In order to understand quantitatively the contributions of the correlated quadrupole 
operators with various ranges of the tensor correlation functions, in Fig.|71we plot the ratio of the 




(50) 




(51) 
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transition matrix elements for the two-body term of the correlated operator against the one with 
the bare operator for each RPA eigenvalue, i.e. 

(vllQ^IIO) 



(52) 



<v||Q||0) 

For the tensor correlator with the constraint ' =0.09 fm , corrections to the bare operator are 
rather small. On the other hand, in the case of tensor correlator with longer range (/^ ' =0.2 
fm^), we notice more pronounced relative contributions from the two-body terms. However, their 
small absolute values result in negligible corrections of the overall transition spectra. It is interest- 
ing to note, that this observation is in agreement with the study of effective operators in the no-core 
shell model within the 2%Q. model space, where the B(E2) values are very similar for the bare and 

n 

the effective operator which includes the two-body contributions [52]. 

Due to the high computational effort involved in evaluating the non-diagonal two-body matrix 
elements for the multipole operators, we used a smaller oscillator basis here {Nmax = 8). We have 
verified that the effect of two-body terms of the correlated operators on the transition spectra is 
very small, regardless of the size of the basis. 



IV. MULTIPOLE EXCITATIONS IN THE UCOM-RPA FRAMEWORK 



Among collective modes of excitations, giant resonances have been a very active topic of nu 



29, 



30, 



31 



3M and 



clear physics in the past few decades, both theoretically 1271 

experimentally Is^ dtI . Collective modes are of particular importance for theoretical models, 
because their underlying dynamics provide direct information about different effective nuclear 
interactions and methods employed to solve the nuclear many-body problem. In the following, us- 
ing the UCOM-RPA scheme, we evaluate the excitation energies and transition strengths of giant 
multipole resonances, and probe their sensitivity on the properties of the correlation functions. 



A. Giant monopole resonances 

The isoscalar giant monopole resonance (ISGMR) is a spherically symmetric oscillation or 
compression of the nucleus, i.e. a breathing mode where neutrons and protons move in phase. The 
ISGMR excitation energy is related to the compressibility of nuclear matter Kn,„, which defines 
basic properties of nuclei, supernovae explosions, neutron stars, and heavy ion collisions Issl . 
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Therefore, it is rather important to provide information about ISGMR from various models based 
on different effective interactions. 

The UCOM-RPA calculated ISGMR strength distributions are displayed in Fig.[5]for a series 
of closed-shell nuclei from ^^O to ^"^Pb. Both the HF single-nucleon basis and the RPA transition 
strength are calculated with the correlated Argonne VI 8 interaction, using the standard ranges 
of the central and tensor correlation functions. For comparison, the unperturbed HF response in 
the 0"^ channel is also shown. One can notice that the unperturbed spectra are widely spread in 
the energy region «20-70 MeV, as a direct consequence of the relatively low level density of the 
UCOM-HF single-particle spectra (Figs. [T] and |2l). However, when we include the RPA residual 
interaction, which is attractive in the isoscalar channel, most of the unperturbed strength is pushed 
to lower energies, generating the collective ISGMR mode. For lighter nuclei '^O, "^"Ca, and "^^Ca, 
the ISGMR is fragmented into two or three peaks, whereas for ^°Zr, '^^Sn, and ^ospj, ^j^g ISGMR 
mode is strongly collective, resulting essentially in a single peak. The calculated ISGMR strength 
distributions shown in Fig. [8] are also compared with the available experimental data from (a, a) 
10.00] and (^He,^He) scattering iSl . and with results from the nonrelativistic RPA [63] and 
relativistic RPA 11641] based on a new interaction with density-dependent meson-nucleon couplings 
(DD-ME2) |0]. The UCOM-RPA excitation energies are in general very close to the values from 
the other studies. In nuclei where the breathing mode is well established ( ^"Zr, ^"*^Pb), the centroid 
energies for the standard correlator are slightly overestimated. The small discrepancies (a: 1-3 
MeV) obtained for the standard correlator could originate from the missing long-range correlations 
beyond the simple mean-field level, the missing genuine three-body interaction, and the two-body 
approximation in the UCOM method, as discussed in Sec. Ill Al 

Next, we address the question to which extent the UCOM-RPA transition spectra are sensi- 
tive to the ranges of the correlators used in the unitary transformation to construct the correlated 
interaction. In Fig. |3 the calculated ISGMR strength distributions are displayed for the corre- 
lated Argonne VIS interaction with different constraints on the range of the tensor correlator, 
' =0.07, 0.08, and 0.09 fm^. One can observe that decreasing of this range systematically 
pushes the transition strength towards lower energies. In particular, by decreasing the range of the 
tensor correlator, i.e. its constraint from /^^~^'^~°^=0.09 fm^ towards 0.07 fm^, the excitation en- 
ergy of ISGMR lowers by «4 MeV. By lowering the range of the tensor correlation functions, the 
density of single-particle spectra increases and the agreement of the ISGMR excitation energies 
with experimental data is improved. 
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B. Giant dipole resonances 



The isovector giant dipole resonances (IVGDR) have recently been studied very extensively, in 
parallel with the renewed interest in the low-lying dipole strength in nuclei away from the valley 



of stability ||66 , 



67, 



68, 



30, HQ- 



Here we employ the UCOM-RPA to evaluate the IVGDR 
strength distributions in light nuclei ^^O, '^''Ca, and "^^^Ca, using the correlated Argonne VI 8 inter- 
action with different constraints on the ranges of the tensor part of the correlator, ' =0.07, 
0.08, and 0.09 fm^ (Fig.fTOb. In lighter nuclear systems, UCOM-RR\ provides a collective char- 
acter of IVGDR, which is distributed over several dominant peaks. We have also extended our 
study of IVGDR to heavier nuclear systems ^"Zr, |^^Sn, and 2°^Pb (Fig.[TT]). The calculated dipole 



response is compared with experimental data ||71 



73 



74 



7511 and with the theoretical excitation 



energies from the relativistic RPA ll64ll based on DD-ME2 interaction [65]. In all nuclei un- 
der consideration, the resulting IVGDR strength distributions display rather wide resonance-like 
structures. The decrease in the range of the tensor correlator, i.e. its constraint /^^~''^~°^=0.09 fm^ 
towards 0.07 fm^, results in lower IVGDR peak energies by ^2-3 MeV 

The calculated IVGDR strength distributions systematically result in higher excitation energies 
than the values from other studies. For the short- range correlator (4^-^'^-°^=0.07 fm^), the cal- 
culated transition strength appears to be in fair agreement with experimental data only for ^^O. 
However, for other nuclear systems, in comparison with experimental data and other theoretical 
studies, the UCOM-RPA overestimates the IVGDR centroid energies by ^^3-7 MeV. This differ- 
ence can serve as a direct measure of the missing correlations and three-body contributions in the 
UCOM-RPA scheme. Inclusion of the three-body interaction and long-range correlations beyond 
the simple RPA method, would probably to a large extent resolve this discrepancies with the other 
studies. Therefore, the 1" channel is particularly convenient to probe the effects of the missing 
correlations and three-body contributions. 



C. Giant quadrupole resonances 

Giant quadrupole resonances are comprised of 2fito ph configurations coupled by the resid- 



ual interaction 11561] . In addition to the resonance-like structure corresponding to the isoscalar 
giant quadrupole resonance (ISGQR), in the isoscalar channel the interaction also generates pro- 
nounced Ofia> low-lying quadrupole states, which have been a subject of various recent theoretical 
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studies ||37l|7^|77l|78|,LZ9|]. The low-lying 2"^ states provide valuable information about the prop- 
erties of the effective interaction 0,0]. 

The isoscalar and isovector quadrupole transition strength distributions are displayed in Fig.fT^ 
for representative cases of medium-mass C^^Ca) and heavy nuclei (^°^Pb). The UCOM-RPA cal- 
culations are based on the correlated Argonne VI 8 interaction with the standard ranges of central 
and tensor correlation functions. The unperturbed HF response is also shown, representing a broad 
distribution of the strength in the region «2-90 MeV which directly reflects the properties of the 



UCOM-HF single particle-spectra (Sec. IIIDI) . When the residual interaction is included in RPA, 
the correlated realistic NN potential generates pronounced low-lying quadrupole excitations in the 
isoscalar channel, and strongly collective states at higher energies corresponding to the ISGQR. 
For "^^Ca, the ISGQR is fragmented into two main peaks in the region «22-30 MeV, whereas for 
the case of ^°^Pb, a single highly collective peak is generated at 20. 1 MeV. In comparison to the 
unperturbed spectra, the transition strength in the isoscalar channel is systematically pushed to- 
wards lower energies. On the other hand, in the isovector channel which is repulsive, the strength 
distributions are in general moved towards the energies above the unperturbed response, resulting 
in broad structures of the transition strength. 

In Fig. [121 we show the UCOM-RPA isoscalar quadrupole transition strength distributions for 
"^"Ca, ^"Zr, and ^"*^Pb. The correlated Argonne VI 8 interaction is employed, with various con- 
straints on the range of the tensor correlator, /^^-^•^-°^=0.07, 0.08, and 0.09 fm^. The calculated 
ISGQR strength distributions are also compared with the experimental da ta |82l] . and the rela- 
tivistic RPA calculations I64I1 with density-dependent DD-ME2 interaction J^]. In all cases, the 
residual interaction constructed from the correlated realistic NN interaction is attractive in the 
isoscalar channel, generating strongly collective peaks corresponding to ISGQR. In addition, in 
the case of ^°Zr, and ^°^Pb, UCOM-RPA also results in the pronounced low-lying quadrupole 
states. The energy of the low-lying quadrupole state in ^°Zr is slightly higher than the experimen- 

I n ' * 

tal value (2.18 MeV Jsj), and for ^'^pb it is rather weU described (4.08 MeV [84]). However, 
the correlated realistic NN interaction is not sufficient for a quantitative description of the ISGQR 
excitation energy. Even the short ranged tensor correlator (/^ ~ ' ~ =0.07 fm ), overestimates 
the experimental values by approx. 8 MeV. By decreasing the range of the tensor correlator, the 
quadrupole response is systematically pushed towards lower energies. In comparison with the 
IVGDR, the quadrupole response has a similar dependence on the range of the tensor correlator. 
For '^''Ca and the tensor correlator ranges determined by 4^"'''^"*'^ =0.07, 0.08, and 0.09 fm^, the 
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ISGQR centroid energies read 25.1, 26.2, and 27.1 MeV, respectively. In the cases of heavier 
nuclei, these differences are smaller: for ^"^Pb, the centroid energy lowers by 1 .2 MeV when go- 
ing from the correlator with ll^~^'^~^^=0.09 inv" towards /^^~''^~*'^=0.07 fm^. From a comparative 
study of low-lying excitations and giant resonances, we can test the sensitivity of different parts of 
the quadrupole transition spectra on the ranges of the tensor correlator (Fig. [131. Whereas the low- 
energy Ofia> excitations depend only weakly on " ' " , the ISGQR appears to be more sensitive 
on variations of the tensor correlator range. 

From our results on the ISGMR, IVGDR and ISGQR nuclear response we can conclude on 
the properties of the correlated NN interactions, used here as effective interactions in the RPA 
calculations. The agreement achieved between the calculated and experimental properties of the 
ISGMR indicates that the correlated NN interaction corresponds to realistic values of the nuclear 
matter (NM) incompressibility. It has been demonstrated in the past that, within relativistic and 
non-relativistic RPA, the energies of the dipole and quadrupole resonances, on one hand, and 
the value of the effective mass corresponding to the effective interaction used, on the other, are 

nn 

correlated II85L 18611 . In particular, the relativistic RPA without density-dependent interaction terms, 
based on the ground state with a small effective mass and relatively high compression modulus, 
resulted in systematically overestimated energies of giant resonances MM- The discrepancies 
between UCOM-RPA calculations and experimental data for multipole giant resonances, as well 
as the low density of single-nucleon UCOM-HF states, suggest that the respective effective mass 
is too small. These observations are consistent with an exploratory UCOM-HF calculation in NM. 



V. SUMMARY 



In the present study, the correlated interaction Vucom based on the Argonne VIS potential 
is employed in fully self-consistent RPA calculations across the nuclear chart. The short-range 
central and tensor correlations induced by the NN interaction are treated within the Unitary Corre- 
lation Operator Method. The precision and self-consistency of the present UCOM-RPA approach 
are tested in the cases of separation of the spurious center-of-mass motion, and in recovering 
the sum rules. It is illustrated that this method provides a highly accurate separation of spurious 
components from the physical spectra, and the sum rules are exhausted up to +3%. A consistent 
implementation of correlated transition multipole operators, even for long-range tensor correlators, 
results in similar response as the bare operators. 
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The Hartree-Fock plus UCOM-RPA framework is employed in studies of giant resonances both 
in light and heavier closed-shell nuclei. The sensitivity of global properties of giant resonances on 
the range parameters of the tensor correlators are systematically studied. The excitation energies of 
giant resonances slowly decrease as the range of the tensor correlator is reduced from the standard 
value ='''^=°^=0.09 fm^ towards 0.07 fm^ 

In comparison to the experimental data, the UCOM-RPA scheme with the standard correlator 
results in slightly higher resonance energy of the breathing mode (~ 1-3 MeV), whereas it overesti- 
mates the excitation energies of IVGDR and ISGQR by « 3-8 MeV. Decreasing of the range of the 
tensor correlator which improves the description of ISGMR is not sufficient to reproduce experi- 
mental data on excitations with higher multipolarities (IVGDRJSGQR). The increased excitation 
energies of giant resonances are related to the rather wide HF single particle spectra which are 
used as a basis for the RPA configuration space, i.e. the small value of the effective mass. Tensor 
correlations with shorter range increase the single-particle level density, improve the description 
of nuclear radii, and result in a systematic shift of the giant resonances towards lower energies. 

The correlated realistic interactions are sufficient to generate collective excitation modes, 
but for an accurate description of experimental data on peak excitation energies and transition 
strengths, the UCOM-RPA approach should be further extended. In particular this would include 
the long-range correlations beyond the simple RPA model, i.e., the inclusion of more complex 
configurations, e.g. the second RPA with coupling of ph with Iplh configurations. The present 
scheme could also be improved by including the self-consistent coupling to two-particle-hole 
phonons within the framework of the dressed RPA IstI . This approach is advantageous since 
it accounts for both, the effects of nuclear fragmentation and the RPA correlations due to two- 
phonon fluctuations in the ground state. 

Apart from long-range correlations, three-nucleon forces might play an important role for the 
quantitative description of single-particle spectra and collective response. Here, we have used 
correlated two-nucleon interactions which are able to describe binding energies in no-core shell- 
model or many-body perturbation theory without any supplementary three nucleon force. This is 
possible due to a cancellation between the attractive genuine three-nucleon force, which accom- 
panies the bare Argonne VI 8 potential, and the repulsive induced three-body force resulting from 
the unitary transformation of the Hamiltonian. This cancellation works surprisingly well for the 
binding energies over a large mass range. However, this does not hold for other observables, e.g. 
the charge radii. Preliminary results using a phenomenological three-nucleon force in the Hartree- 
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Fock framework indicate that radii and single-particle spectra can be improved significantly. As a 
consequence an improvement of the calculated excitation energies of giant resonances with respect 
to the experimental data can be expected. 

The inclusion of the missing correlations and residual three-body interactions in the UCOM- 
RPA scheme, that should improve the description of low-lying excitations and giant resonances on 
a quantitative level, will be the subject of our future work. In addition, the UCOM-RPA correla- 
tions due to low-lying states and collective excitations, may also provide the long-range correla- 
tions which are necessary to improve the HF ground state. 
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FIG. 1: The UCOM-HF neutron and proton single particle spectrum for ^^O, along with the corresponding 



HP spectra based on the low-momentum NN potential Viow-k Il25ll . phenomenological HF with the SIII 



Skyrme-type interaction 



111 



, relativistic mean field theory with the NL3 effective interaction 



and the 



experimental spectra ll46ll . The UCOM-HF calculations are based on the correlated Argonne V 1 8 interaction 



with the standard constraints on the ranges of the central and tensor correlators (Sec. Ill At . 
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FIG. 2: The same as Fig.lH but for ^^Ca. 
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FIG. 3: The number of ph configurations in the UCOM-RPA (right scale) and the energy of the spurious 
center-of-mass state (left scale) as a function of the cut-off for the maximal energy of the ph transitions. 
Two sample nuclei are considered, '^O and ^°*^Pb. The correlated Argonne V18 interaction is employed, 
using the standard correlators (Sec. Ill Ab . 
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FIG. 4: The UCOM-RPA isoscalar dipole transition spectta for '^'^Ca, by employing the isoscalar dipole 
operator without (thin solid Une) and with (wide grey line) the correction term for the spurious center-of- 
mass motion (Argonne V18, /i^=^'^"°^=0.09 fm^). 
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FIG. 5: The cumulative sum of the UCOM-RPA energy weighted strength for the isoscalar monopole 
(upper panel) and the quadrupole response (lower panel), in comparison with the energy weighted sum 
rules (horizontal dotted lines) (EHli,<E3, for '^O, "^^^Ca, ''^Zr, and ^o^pb (Argonne V18, //^''^^''^-O.OQ fm^). 
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FIG. 6: The isoscalar monopole (upper panel) and quadrupole (lower panel) strength distributions for °0, 
calculated with the bare and correlated transition operators. The correlated Argoime V18 interaction is 

(S=l 7^=0) 

used, with the constraints ' -0.09 and 0.2 fm^ for the range of the tensor correlator (Nmax - 8 for 
the harmonic oscillator basis). 
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FIG. 7: The absolute value of the ratio of the transition matrix elements for the two-body terms of correlated 
operator against those of the bare operator. The tensor correlators are constrained by ' -0.09 fm^ 
and //=^'^=°^=0.2 fm^ (Argonne V18, N,„ax = 8). 
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FIG. 8: The UCOM-RPA transition strength distributions of the ISGMR (solid line), in comparison with the 
unperturbed HF spectra (dashed Une). The correlated Argonne V18 interaction is used, with the standard 



constraint (/l^ ''^ "^=0.09 fm^) on the range of the tensor correlators. The experimental data 5S , 6C , 



6211 . results from the nonrelativistic (Drozdz et 



al.) El, 



and the relativistic RPA (DD-ME2) M , 



u 

65ll are 



denoted by arrows. 
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FIG. 9: The calculated UCOM-RPA strength distributions of the ISGMR (Argonne VI 8), using different 
restrictions on the range of the tensor correlator (/^^"'•^""^^O.OV, 0.08, and 0.09 fm^). Results from the 
other studies, denoted by arrows, are the same as in Fig. |S] 
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FIG. 10: The UCOM-RPA strength distributions for the IVGDR in ^^O, "^^Ca, and "^^Ca. The calculations 
are based on the correlated Argonne VI 8 interaction, using different constraints on the tensor correlator 



range (/ 



.(S=i,r=o)_ 



0.07, 0.08, and 0.09 fm^). The experimental }1J\ and theoretical IVGDR energies from 



the relativistic RPA (DD-ME2) ||64 , 



6511 are denoted by arrows. 
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FIG. 11: The same as Fig. [13 but for ^^Zni^Zsn, and ^o^Pb. The experimental data in 



relativistic RPA (DD-ME2) energies ll64ll65ll are shown by arrows. 
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FIG. 12: The ISGQR and IVGQR strength distributions for '^'^Ca and ^"'^Pb, in comparison with the unper- 
turbed HF response (Argonne V18, if =^''^^^^^0.09 fm\ 
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FIG. 13: The UCOM-RPA strength distributions for the ISGQR in ^OCa, ''"Zr, and ^o^pb. The correlated Ar- 

(5 = 1 T=0) 

gonne V18 interaction is used, with different ranges of the tensor correlator (/^ -0.07, 0.08, and 0.09 



fm^). The experimental ISGQR excitation energies 



RPA energies 164 



El, 



the low-lying states |83 



8411 . and the relativistic 



6511 are denoted by arrows. 
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